Spin fluctuations and superconductivity in a 3D tight-binding model for BaFe2As2 
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Despite the wealth of experimental data on the Fe-pnictide compounds of the irFe2As2-type, K 
— Ba, Ca, or Sr, the main theoretical work based on multiorbital tight-binding models has been 
restricted so far to the study of the related 1111 compounds. This can be ascribed to the more 
three dimensional electronic structure found by ab initio calculations for the 122 materials, making 
this system less amenable to model development. In addition, the more complicated Brillouin zone 
(BZ) of the body-centered tetragonal symmetry does not allow a straightforward unfolding of the 
electronic band structure into an effective IFe/unit cell BZ. Here we present an effective 5-orbital 
tight-binding fit of the full DFT band structure for BaFe2As2 including the kz dispersions. We 
compare the 5-orbital spin fluctuation model to one previously studied for LaOFeAs and calculate 
the RPA enhanced susceptibility. Using the fluctuation exchange approximation to determine the 
leading pairing instability, we then examine the differences between a strictly two dimensional model 
calculation over a single kz cut of the BZ and a completely three dimensional approach. We find 
pairing states quite similar to the 1111 materials, with generic quasi- isotropic pairing on the hole 
sheets and nodal states on the electron sheets at /c^ = which however are gapped as the system is 
hole doped. On the other hand, a substantial kz dependence of the order parameter remains, with 
most of the pairing strength deriving from processes near kz = tv. These states exhibit a tendency 
for an enhanced anisotropy on the hole sheets and a reduced anisotropy on the electron sheets near 
the top of the BZ. 

PACS numbers: 74.70.Xa,74.20.Pq,74.20.Rp 



INTRODUCTION 



The discovery of superconductivity in the Fe-based 
pnictide and chalcogenide compounds has stimulated a 
tremendous research effort in many areas of solid state 
physics and chemistry. The materials initially discov- 
ered by the Hosono group, LaOFeP^ and fluorine doped 
LaOFeAs^, belong to a class of iron pnictides that are 
commonly referred to as the 1111 structures. More re- 
cently, a great deal of attention has been devoted to 
the preparation and investigation of materials where the 
FeAs layers are separated by a single cation only, the 
so-called 122 structures^. Here both hole-doping by re- 
placing Ba in part by K as well as electron-doping due 
to a fractional substitution of Fe by Co have proven suc- 
cessful in suppressing the spin-density wave (SDW) for- 
mation in favor of a superconducting ground state^i^. Al- 
though the maximum critical temperature of Tc = 38K 
in the 122-systems is smaller than in the related 1111- 
materials, the possibility of growing high quality sin- 
gle crystals with relatively good surfaces make them op- 
timal candidates for comprehensive experimental stud- 
ies. Angle-resolved photoemission electron spectroscopy 
(ARPES) measurements performed on high quality single 
crystals of Bao.6Ko.4Fe2As2r'— have been very influen- 
tial, revealing the position, shape and size of the Fermi 
surface pockets that are in qualitative agreement with 
band structure calculations^iii. In addition, the ARPES 



experiments also claim to resolve the size and momentum 
space distribution of the superconducting gap, showing 
at least two distinct values of the order parameter and 
a nearly isotropic gap size distribution along the indi- 
vidual Fermi surface sheets. These observations can in 
principle be understood in terms of the formation of a 
sign-changing s-wave state generated by the exchange of 
spin fluctuation^ii^. Such a gap structure is also sup- 
ported by neutron scattering experiments on these com- 
pounds, which find a resonance emerging in the super- 
conducting state at a wave vector that corresponds to 
the separation between hole and electron pocket^"—. 
Despite these promising results, the symmetry of the su- 
perconducting order parameter in these materials is still 
controversial, and many experiments imply the existence 
of low-energy quasiparticle excitations, denoting the pos- 
sible existence of nodes. These include NMR^^^^, super- 
fluid densit}i2^ — , thermal conductivit y^^ and Raman 
light scattering^. The existence of nodes is also sug- 
gested by multiorbital calculations for the 1111 materi- 
g^lg34,35^ although a transition from such a nodal state to 
a fully gapped state has also been described in such a 
framework^—. While these works indicate a sensitivity 
of the superconducting state to details of the electronic 
structure, particularly the position of the pnictogen, they 
have not been discussed in the framework of the 122 ma- 
terials. 

In particular, it is important to ask if the greater three- 
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FIG. 1. (Color online) Sketch of the Brillouin zone of the 
lA/mmm crystal symmetry (a) and of the large effective BZ 
corresponding to the IFe/unit cell (b). The blue line shows 
the two paths in the IFe/unit cell BZ that have to be folded 
by the reciprocal lattice vector T = (tt, tt, tt) (red arrow) to 
give the corresponding path in the 2Fe/unit cell BZ of the 
P4 / nmm symmetry. 

dimensionality of the electronic structure in these mate- 
rials has qualitative effects on the pair state. Certainly 
the effects of near nesting of the hole and electron sheets 
of the Fermi surface, which are said to be responsible for 
the stabilization of the s± state, must be expected to 
change for those sheets which are strongly /c^ -dispersive. 
One important aspect in this regard is the role of doping 
since different impurities may not only provide carriers 
for the FeAs planes but can also alter the hybridization of 
states between layers, as argued for the case of Co doping 
by Kemper et al^. In general, the more three dimen- 
sional character and the resulting coupling between the 
individual FeAs layers in the 122 compounds requires us 
to revisit the existing models and to reassess the univer- 
sality of the results obtained so far. 

II. TIGHT-BINDING FIT OF THE LDA BAND 
STRUCTURE 

Compared to the other pnictide and chalcogenide su- 
perconductors as e.g. the 1111 materials (LaOFeAs), 
the 111 systems (LiFeAs), or the binary 11 compounds 
(FeSe), there is no straightforward way to describe the 
electronic structure of the 122 materials (BaFe2As2) in 
an effective 5-orbital Fe model. This has several reasons: 
first of all the Ba atoms, forming a spacing layer between 
the FeAs planes, contribute significantly to the inter layer 
hopping. Secondly, we have only one Ba atom per unit 
cell, but two Fe and two As atoms, therefore a descrip- 
tion in an effective model based on a IFe/unit cell BZ 
can only be successful if the Ba bands are integrated out 
from the beginning. Finally the basic unit cell in the 122 
materials is not a simple tetragonal unit cell but a body- 
centered tetragonal unit cell, with non-Cartesian recipro- 
cal lattice vectors and a complicated Brillouin zone (see 
Fig.d]). Despite these difficulties it is in principle possible 
to define an effective band structure in the IFe/unit cell 
BZ starting from a purely Fe-based fit of the full DFT 
band structure^. 

We have calculated the band structure for the 
BaFe2As2 parent compound making use of the density 
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FIG. 2. (Color online) The paramagnetic DFT band struc- 
ture (full line) and a Wannier fit (crosses) of the 10 bands 
in the vicinity of the Fermi surface onto the Fe-Sd orbitals 
(a). The 5-orbital tight-binding fit (colored points) of the 10- 
orbital Wannier fit (black points) with a color coding of the 
main orbital contributions (b). The colors correspond to dxz 
(red), dyz (green), dxy (blue), dx2_y2 (orange), and 6/3^2 
(magenta) . 



functional theory (DFT) in a plane wave basis set with 
ultrasoft pseudopotentials as provided in the Quantum 
ESPRESSO packageM. Here we have used the lattice 
constants as well as the internal coordinates tabulated 
in Ref. |4a with a = 3.9625 A, c = 13.0168 A, and 
ZAs = 0.3545. 

The calculations were performed for the body-centered 
tetragonal unit cell corresponding to the R/mmm sym- 
metry of the crystal but we have plotted the bands along 
the high symmetry lines of a corresponding simple tetrag- 
onal unit cell to facilitate the comparison with the band 
structure of the 1111 materials. In the next step, we pro- 
jected the bands in the vicinity of the Fermi energy on the 
Fe-3d orbitals using maximally localized Wannier func- 
tions (MLWF) following the method of Marzari and Van- 
derbilt^. The bands were disentangled by minimizing 
the spread of the Wannier functions. Except for a band 
with mainly Ba character that approaches the Fermi en- 
ergy between F and X this projection reproduces the full 
DFT band structure very accurately as can be seen from 
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FIG. 3. (Color online) The partial density of states of the 
5-orbital tight-binding fit, using the same color coding as in 
Fig.Elb. 

Fig. [2] a. Finally, we fitted the Wannier bands with a 
5-orbital tight-binding Hamiltonian, unfolding the small 
2Fe/unit cell BZ to a large IFe/unit cell effective BZ. The 
Hamiltonian is given as 

Ho = Y.Yl (^^"(1^) + «-^-") ^L(k)rf„a(k) (1) 

kcr rnn 

where dj^ ^ (k) creates a particle with momentum k and 
spin a in the orbital m. The kinetic energy terms 
£,mnO^) together with the parameters of the 5-orbital 
tight-binding fit are listed in the appendix. The orbital 
resolved density of states shown in Fig. [3] reveals dom- 
inant contributions to the total density of states at the 
Fermi level from the Fe dxz/dyz and dxy orbit als. In 
Fig. H] we show the Fermi surface pockets at two different 
kz cuts of the BZ, where the colors encode the main or- 
bital contributions to the respective band. As discussed 
in Kemper et al.^^^ the important orbital matrix elements 
in the pairing interaction enter as |a^(k)p, where t de- 
notes the orbital and v the band index. Plots of |a^(k)p 
for kz = and kz = tt are shown in Fig. [5] a and b, re- 
spectively. This figure is similar to Fig. HI but contains 
more detailed information, since it shows all of the orbital 
contributions, while Fig. |4] shows only the largest orbital 
contribution on a given part of each Fermi surface. Here 
we note that the orbital composition of the Fermi sur- 
faces at kz = is similar to the one found for the 1111 
materials (e.g. compare with Fig. 5b in Ref. |35[ ), while 
at kz = 7T the orbital composition changes substantially. 
Here the inner hole pocket around F (labeled as ai) is of 
predominantly dxy character while the outer hole pocket 
around F (labeled as 0^2) has in addition to the dxz/yz 
contributions significant involvement of the dx'^-y^ or- 
bit als. From kz = to kz = 7T one also finds a change 
in the orientation of the ellipticity of the electron pock- 
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FIG. 4. (Color online) The main orbital contributions to the 
Fermi surfaces of the hole doped compound with (n) = 5.9 at 
kz = (a) and kz = tt (b) using the same color coding as in 
Fig.[2]b. 



ets Pi and P2 without a qualitative change of the orbital 
composition. Recently, a similar procedure was used by 
Miyake et al^ to determine an effective tight-binding 
model ^mn(k) for BaFe2As2 based on a DFT calculation. 
Their results appear to agree roughly with ours, with the 
exception that they used a different choice for the phase 
of the orbital basis, leading to a discrepancy in the sign 
of some of the dispersions. This difference should not 
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FIG. 5. (Color online) The orbital composition of the Fermi 
surface sheets for a hole doped compound ((n) = 5.9) as given 
by I at (k) 1 2 for kz = (a) and kz = tv (b). The orbital 
contributions are shown as a function of the winding angle 
a starting with the rightmost point on each Fermi surface 
sheet. 



affect the eigenenergies, however. 



III. 3D MULTIORBITAL SUSCEPTIBILITY 

In this section we discuss the differences between a sus- 
ceptibility calculated for a fixed value of kz neglecting the 




FIG. 6. (Color online) Fermi surface mesh of the hole doped 
compound applied to the calculation of the pairing functions. 
Here we used 24 x 10 k-points for every Fermi surface sheet 
with ai (red), a2 (blue), /32 (green), and 7 (yellow). 



kz dispersion of the energy bands and a complete 3D cal- 
culation of the susceptibility, taking the full momentum 
dependence of the band structure into account. Using 
the notation of Ref. |47| we write the non-interacting sus- 
ceptibilities as 

\<.,iuJn 

(2) 

where TV is the number of Fe lattice sites, /3 = 1/T is 
the inverse temperature, ujn are the fermionic and oOm 
the bosonic Matsubara frequencies in the imaginary time 
formalism, and s, t, p, and q are indices denoting the Fe- 
3d orbit als. For the full 3D susceptibility the momentum 
sum runs over kx^ ky^ and /c^, while for the 2D calcula- 
tions kz is kept fixed and the susceptibility is only evalu- 
ated dit Qz = 0. For the integration we use a 64 x 64 x 20 
/c-mesh and we notice only negligible finite size effects, 
that show up as weak oscillations at small q. The spec- 
tral representation of the Green's function is given as 



Gsp{k,iuj„) 



E 

/i 



a^(k)ar(k) 



(3) 



where the matrix elements a^(k) = (5|/ik), connecting 
the orbital and the band space, are determined by a di- 
agonalization of the intra- and interorbital dispersions of 
the tight-binding Hamiltonian given in Eq. [H Now we 
calculate the retarded susceptibility as 

a^(k)a^*(k)<(k + q)ar(k + q) 



N 

x[/(^.(k- 



■E,{k + (i)-E^{k) + iO+ 



■q))-/(^^(k))] 



(4) 



It is evident that for a system without kz dispersion the 
susceptibility does not depend on Qz and the sum over kz 
can be neglected. Finally, we take the orbital dependent 
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FIG. 7. (Color online) 2D susceptibility: The real part of the 
RPA enhanced susceptibility Xrpa(^) as a function of the in- 
plane momentum transfer calculated in two dimensions for a 
single value of /c^, /c^ = (a,c,e) and kz = tv (b,d,f) for a 
hole doped compound with (n) = 5.9. For (a) and (b) we 
have used U — 0.65 and J = 0, while for (c) and (d) we have 
used U = 0.55 and J = 0.25U. In panels (e) and (f) the 
susceptibility is shown along the main symmetry lines with 
U = 0.65, J=0 (red), and U = 0.55, J = 0.25U (blue). 



FIG. 8. (Color online) 3D susceptibility: The real part of the 
RPA enhanced susceptibility xrpa(^) as a function of the in- 
plane momentum transfer for two different values of qz, Qz = 
(a,c,e) and qz — (b,d,f) for a hole doped compound with 
(n) = 5.9. For (a) and (b) we have used U = 1.1 and J = 0, 
while for (c) and (d) we have used U — 0.8 and J — 0.25[/. 
In panels (e) and (f) we use the same coloring scheme as in 
Fig. El 



interactions into account by defining the RPA enhanced 
spin susceptibility as 

(Xl )stpq = Xstpq (5) 

where we sum over repeated indices. Here the interaction 
matrix is nonzero only for 

where a ^ h and the definitions of the intraorbital repul- 
sion /7, interorbital interaction U' ^ Hund's rule coupling 
J, and pair hopping energy J' in terms of a general inter- 
action Hamiltonian are given in Ref. 113 and are related 
to the notationin Ref. U = U, U' = V ^ J/ 4., 

J = J/2, and J' = J^ In the spin rotational invariant 
case the interaction parameters are connected by J = J' 
and U' -\- J' = U — J. The charge susceptibility 

iXo'^^) stpq can be derived similarly with a different in- 
teraction matrix with the components 

^aaaa = ^5 ^bbaa = "^U' — J, t/^6a6 = 2J—U', U^^j^^ = J' . 



In the following we will distinguish between sets of inter- 
action parameters with and without finite Hund's rule 
coupling (and corresponding pair hopping) and we will 
adjust the interaction parameters to be close to the su- 
perconducting instability. Note that the ratio oi U / J 
taken here for J > cases is similar to that found by 
Miyake et al^ from ah initio calculations, but that the 
overall scale is smaller. This renormalized RPA interac- 
tion scale is familiar from one-band interacting models in 
the cuprates^^. In Fig. [71 we show the real part of the 
RPA enhanced spin susceptibility 

X^(^) = \ Y^{Xi)sspp{^. 0) (6) 

sp 

calculated over a 2D Fermi surface corresponding to a cut 
through the 3D Fermi surface at fixed kz=^ or tt, and for 
two different sets of interaction parameters. Here we find 
that the susceptibility calculated at /c^ = shows only a 
small incommensurate enhancement of the main scatter- 
ing peaks at q = (tt, 0), while the susceptibility calculated 
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FIG. 9. (Color online) 2D susceptibility: The real part of the 
RPA enhanced susceptibility Xrpa(^) as a function of the in- 
plane momentum transfer calculated in two dimensions for a 
single value of A:^, kz = (a,c,e) and kz = tv (b,d,f) for an 
undoped compound with (n) = 6. For (a) and (b) we have 
used U = 0.7 and J = 0, while for (c) and (d) we have used 
U = 0.6 and J = 0.25U. In panels (e) and (f) we use the 
same coloring scheme as in Fig. [T] 



FIG. 10. (Color online) 3D susceptibility: The real part of 
the RPA enhanced susceptibility xrpa(^) as a function of the 
in-plane momentum transfer for two different values of qz, 
Qz — (a,c,e) and qz — (b,d,f) for an undoped compound 
with (n) =6. For (a) and (b) we have used [7=1.1 and J = 0, 
while for (c) and (d) we have used U — 0.9 and J — 0.25[/. 
In panels (e) and (f) we use the same coloring scheme as in 
Fig. [3 



at /c^ = TT is commensurate and strongly peaked for the 
same value of the momentum transfer. This enhancement 
at /c^ = TT can be explained by the additional intraorbital 
scattering channel between the electron and hole pockets 
due to a major dxy contribution on the ai FS sheet at 
kz = TT that is absent at /c^ = (see Fig. 2]). Compar- 
ing the susceptibility along the main symmetry lines in 
panel (e) and (f), we find that for the kz = tt cut the 
non-zero Hund's rule coupling J = 0.25/7 together with 
the related pair hopping J' = J leads to a strong en- 
hancement of the scattering peak at q = (tt, 0) exceeding 
the scattering peak at q = (7r,7r). By contrast, an exclu- 
sive increase of the intra-orbital interaction U without 
Hund's rule coupling J = 0, chosen to produce a compa- 
rable enhancement of the scattering peak at q = (7r,0), 
also enhances the scattering peak at q = (7r,7r). 

Comparing the results for the two different kz values 
we therefore expect that for the full 3D calculation, con- 
tributions from kz = TT will dominate the total suscepti- 
bility but the structures will be less sharp as a result of 
the kz averaging process. In Fig. [51 we have now studied 
the full 3D susceptibility including the explicit Qz depen- 



dence and summing over kz. As expected from our two 
dimensional study the main peak structure is similar to 
the dominating kz = tt cut of the susceptibility. But we 
note that for J = the well separated scattering peak at 
q = (tt, tt) in the kz = tt susceptibility is already broad- 
ened to a plateau-like structure with its main weight 
shifted to an incommensurate position around q = (tt, tt). 
For J > the substructures in the scattering peaks van- 
ish and the susceptibility shows broad but well developed 
(tt, 0) scattering peaks at a commensurate position. Note 
this result is quite similar to the commensurate normal 
state neutron scattering intensity observed by Inosov et 
ali^, which was difficult to understand in the 2D calcu- 
lations for the 1111 Fermi surface in Ref. |35|. Here we 
find that the peaks in Re x(Q7^ = 0) tend to correspond 
to those of Im x{^^^) small cj, suggesting that the 
"averaging" of the susceptibility due to the 3D disper- 
sion in the 122 materials is sufficient to account for the 
commensurate response. 

In Figs. [9] and [TOl we show the equivalent results to 
Figs. [3 and [8] but for zero doping. Compared to the hole 
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FIG. 11. (Color online) 2D pairing functions, J = 0: The 
leading (upper row) and subleading (lower row) pairing func- 
tion for the hole doped compound ((n) = 5.9) plotted along 
the Fermi surfaces at two different kz cuts. The pairing func- 
tions are shown in the order ai, a2, 7, and /32 running 
counter-clockwise around each Fermi surface sheet with the 
rightmost point as the starting point on each sheet, except the 
(32 pocket where the plots start with the uppermost point. 
The calculations were performed for U = 0.65 and J = 
and the eigenvalues are A = 0.022 (s-wave) and A = 0.015 
(d-wave) for kz = 0, and A = 1.615 (s-wave) and A = 0.377 
((i-wave) for kz = tv. 



doped case, we see that the peaks at (tt, 0) are consid- 
erably suppressed (see Fig. [9]), and new incommensurate 
peaks along the (tt, 0)-(0, tt) line appear. In the fully 
integrated 3D result, Fig. [TOl by contrast, the nearly 
commensurate response at (tt, 0) is recovered, but is sup- 
pressed relative to an incommensurate ridge of response. 



IV. PAIRING SYMMETRY 

Again following the notation in Ref. 113 we define the 
singlet pairing vertex in the fluctuation exchange approx- 
imation^!^ as 



-/7^X^^(k-k^c^)/7^ 



i/7^X?P^(k-k^c^)/7^ + i/7^ 



(7) 



pstq 



where x?^^ denotes the RPA enhanced spin and Xo^^^ 
the RPA enhanced charge susceptibility. Making use of 
the Kramers-Kronig relation we can proceed further by 
considering only the real part of the static pairing interac- 



FIG. 12. (Color online) 2D pairing functions, J > 0: The 
leading (upper row) and subleading (lower row) pairing func- 
tions for the hole doped compound ((n) = 5.9) plotted as 
before for two different values of kz, calculated for U = 0.55 
and J — 0.25[/. The eigenvalues are A = 0.027 (s-wave) and 
A = 0.015 (d-wave) for kz = 0, and A = 0.638 (s-wave (1)) 
and A = 0.381 (s-wave (2)) for kz = tt. 



tion. This procedure is justified in the regime of weak to 
intermediate correlations considered herein and was dis- 
cussed in detail in a previous work—. If we now confine 
our considerations to the vicinity of the Fermi surfaces 
we can determine the scattering of a Cooper pair from 
the state (k, — k) on the Fermi surface Ci to the state 
(k^, — k') on the Fermi surface Cj from the projected in- 
teraction vertex 

ri,(k,k') = Re^<^(k)at.(-k)ri,p,(k,k',0) 



stpq 

x<f(k')a«': 



(-k') 



(8) 



where the momenta k and k' are restricted to the differ- 
ent Fermi surface sheets with k G and k' G Cj . Defin- 
ing a dimensionless pairing strength functional we calcu- 
late the symmetry function ^^(k) of the leading pairing 
instability from the following eigenvalue problem 



E 



-r^^-(k,k')^c.(k') = XaQaO^) 

(9) 

where '^^(k) = |Vk£^i/(k)| is the Fermi velocity on a given 
Fermi surface. The largest eigenvalue will lead to the 
highest transition temperature and its eigenfunction de- 
termines the symmetry of the gap. For the numerical 
calculation of the hole doped compound we parametrize 
the Fermi surface by a dense mesh of 1200 k values 
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FIG. 13. (Color online) 3D pairing functions, hole doped: 
The leading pairing functions for the hole doped compound 
((n) = 5.9) plotted as before for two different values of kz, 
calculated for [7 = 1.0 and J = (upper row) and [7 = 0.8 
and J — 0.25[/ (lower row). Here the maximum eigenvalues 
are A = 0.956 and A = 1.077, respectively. 
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FIG. 14. (Color online) 3D pairing functions, undoped: The 
leading and subleading pairing functions (d-wave and s-wave) 
for the undoped compound ( (n) — 6) plotted as before for two 
different values of kz, calculated for U — 0.9 and J — 0.25[7. 
Here the eigenvalues are A = 1.14 and A = 0.617, respectively. 



distributed over the 5 different Fermi surface sheets as 
shown in Fig. [6l 

First we study the pairing function at a fixed kz cut of 
the BZ. In order to solve the eigenvalue problem of Eq. [9] 
we first use an effective pairing interaction rij(k, k') cal- 
culated from the 2D susceptibility derived in the previous 
section. Due to the different orbital contributions along 
the Fermi surface sheets at /c^ = and /c^ = tt as well as 
due to the different nesting conditions, reflected in the 
respective 2D susceptibilities, we find for a given set of 
interaction parameters quite different pairing functions 
at the center and on top of the BZ. In Fig. [11] we show 
the pairing functions for /7 = 0.65 and J = for a hole 
doped compound. At /c^ = the rather featureless sus- 
ceptibility (compare Fig. [71 a), without a distinct peak at 
the nesting vector q = (tt, 0), results in a leading 5- wave 
pairing instability with a sign change between the a/7 
and the j3 FS sheets. At /c^ = tt the strong scattering 
peak in the susceptibility at q = (tt, 0) (Fig. [71b) drives 
also an 5-wave pairing state, but the frustration intro- 
duced by the equally strong scattering peak at q = (tt, tt) 
enforces an additional sign change between the a and 7 
sheets and suppresses the gap on the electron pockets. 
Comparing the eigenvalues for the two different kz cuts 
using the same set of interaction parameters we find that 
we are still far from the instability at /c^ = while we 
already have a divergent eigenvalue at kz = tt. 

If we reduce the intra-orbital pairing interaction IJ = 



0.55 and simultaneously add a finite Hund's rule cou- 
pling J = 0.25/7 (Fig. [12]) we find for both kz cuts an s- 
wave symmetry of the leading pairing state without the 
sign change between the a and 7 Fermi surface sheets 
dit kz = TT and a more isotropic gap size along the elec- 
tron pockets. This can be understood in terms of the 
enhanced q = (tt, 0) peak in the susceptibility without a 
simultaneous enhancement of the q = (tt, tt) peak. Here 
the subleading pairing state is either a d-wave {kz = 0) or 
a different 5- wave {kz = tt). Obviously a situation where 
for a fixed set of interaction parameters different pairing 
symmetries at different kz cuts of the Fermi surface might 
be realized does not reflect an energetically favorable so- 
lution: thus it is evident that only a full 3D calculation 
of the pairing state can succeed. However, comparing the 
eigenvalues A for kz = and kz = tt we have confirmed 
that the primary contributions to the pairing come from 
near kz = tt. 

In a next step we study the pairing function calculated 
from a full 3D susceptibility using the complete Fermi 
surface mesh as shown in Fig. [6] Here the susceptibility 
is "averaged" over kz and we have already seen that it 
only weakly depends on Qz. In the upper row of Fig. [131 
we show the leading eigenfunction for U = 1.0 and J = 
dit kz = and kz = tt for a hole doped compound. Here 
the leading pairing state is an extended 5-wave state ex- 
hibiting a higher anisotropy on the electron FS sheets Pi 
and dit kz = than at kz = tt. On the hole pockets 
around F we find at /c^ = a smah but isotropic gap 
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FIG. 15. (Color online) 3D pairing functions, J — 0: The 
leading Alg pairing function (extended s-wave) on the ai, 0^2, 
and 7 Fermi surface sheet for the hole doped compound 
((n) = 5.9). The semi-transparent color mesh visualizes the 
gap on each each of the FS sheets that is also shown without 
the gap for comparison. The calculations were performed for 
U — 1.^ and J = 0. In this figure we have changed the overall 
phase of the gap by -1 from that used in Fig. [13] in order to 
show the nodal structure of the gap more clearly. 



that develops a strong anisotropy towards /c^ = tt, where 
it finally exhibits several sign changes around the a FS 
sheets. We also observe an overall sign change of the gap 
on the hole pockets around F as a function of kz . On the 
7 FS sheet we have a nearly isotropic gap with opposite 
sign compared to the j3 sheets in consequence of the pro- 
nounced (7r,0) scattering peak in the 3D susceptibility. 
For a finite Hund's rule coupling J = 0.25U (lower row 
in Fig. [T3|) we find a very similar pairing state with a 
reduced gap on the a Fermi surface sheets at the top of 
the Brillouin zone. 

In Fig. [151 the gap state from the upper row of Fig. [13] 
is shown along the unwrapped Fermi surfaces for the ai, 
0^2, /^i, and 7 sheets. Here it becomes apparent that 
the sign change of the pairing state as a function of kz 
on the ai Fermi surface leads to nearly horizontal line 
nodes close to the top of the Brillouin zone. On the a2 
sheet, on the other hand, small V-shaped line nodes with 
vertical components are also present. 

For the undoped compound (Fig. [14]) we find a d-wave 
solution as the leading eigenfunction for all parameters 
studied so far. The anisotropic-s solutions are character- 
ized by greater anisotropy on the P sheets with nodes in 
the direction connecting the a and /3 FS sheets. Here we 
notice that the presence of the pair hopping term J > 
again does not change the results qualitatively. This can 



be understood from the susceptibility (Fig. [TQ|) that does 
not show a significant enhancement of the incommensu- 
rate scattering peak near q = (tt, 0) in the presence of 
a finite pair hopping parameter J. Calculations for an 
electron doped compound (not shown here) lead to re- 
sults that are qualitatively similar to the ones for the 
undoped compound with a leading d-wave pairing state 
that is well separated from the subleading sign changing 
5- wave state. 



V. CONCLUSIONS 

In this article we have presented a 5-orbital tight- 
binding fit of the DFT band structure of BaFe2As2 de- 
rived by unfolding the bands of the real BZ into an ef- 
fective IFe/unit cell BZ. Here we found that the orbital 
contributions to the Fermi surface sheets at /c^ = are 
qualitatively different from the ones dit kz = tt where 
the hole pockets around F show a pronounced multior- 
bital composition. Based on this 5-orbital tight-binding 
model we compared the 2D and 3D RPA susceptibilities 
where the latter is calculated by integrating over the full 
Brillouin zone. We showed that the susceptibility is dom- 
inated by contributions from the top and the bottom of 
the BZ and develops well pronounced scattering peaks in 
a certain range of parameters. We also found that the 
Qz dependence of the susceptibility is weak but this does 
not imply that the kz integration of the susceptibility can 
be neglected since it leads to an averaging of the kz de- 
pendent susceptibilities. In particular the commensurate 
nature of the 3D magnetic response may depend on this 
averaging. 

Finally, we studied the pairing functions in the fluc- 
tuation exchange approximation and compared again a 
strictly 2D calculation to a complete 3D calculation. 
Here it becomes obvious that due to the strong kz disper- 
sion and the different orbital composition of the bands in 
the center and on top of the Brillouin zone the pairing 
function changes considerably along the Fermi surface 
cylinders and a 2D description will fail to find the most 
stable pairing state over the full Fermi surface. Within a 
3D approach we showed that for a moderate hole doping 
the existence of the additional hole pocket around (tt, tt) 
favors an extended s-wave state over the d-wave state 
found for the undoped compound. The pairing strength 
for these states is found to arise predominantly from pro- 
cesses with momenta near kz = tt due to the dxy orbital 
character of the hole sheets in this region. The s-wave 
state exhibits a strong anisotropy on the electron pock- 
ets dit kz = that is reduced at kz = tt, while on the 
hole cylinders around F we have a nearly isotropic gap 
dit kz = that develops vertical line nodes and changes 
its sign towards kz = tt. For finite J the gap maximum 
on the hole sheets around F is reduced at kz = tt. Inde- 
pendent of J, the gap on the hole cylinder around (tt, tt) 
is large and isotropic and nearly independent of kz- 

Such states should exhibit responses to external probes 
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quite different than the 2D states which have been dis- 
cussed in the hterature until now. In particular, the 
nodes near the top of the Brillouin zone will contribute 
strongly to low-temperature c-axis transport, and should 
produce, e.g. a strong linear-T term in the penetration 
depth Ac(T) and a weaker behavior in Aa6, as observed 
recently by Martin et al^. 

We emphasize that our primary purpose in this work 
has been to investigate the novel qualitative aspects of 
the spin fluctuation pairing based on the 122-type Fermi 
surface compared to the more familiar 1111-type, in par- 
ticular those aspects driven by the 3D dispersion. As 
such, we do not claim to have fully explored interaction 
parameter space, nor to have chosen those parameters 
appropriate to a particular system, as attempted, e.g. in 
Ref. |43|. Thus while within the limited parameter sets 
we have investigated we have looked at the effects of hole 
doping, it should be understood that these are changes 
in electronic structure which may be controlled by other 
variables which influence the electronic structure, such as 
the pnictogen/chalcogen position, existence of surfaces, 
or presence of disorder. Further work is necessary to see 
which extensions of the current theory are most vital to 
make quantitative comparisons with experiment. 
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Appendix A: Fitting parameters for the 5-orbital 
tight-binding model 

Subsequently we show the dispersions of the 5-orbital 
tight-binding fit of the BaFe2As2 band structure. Here 
we note that unfolding the small 2Fe/unit cell BZ into 
an effective large BZ corresponding to a unit cell with 
only one single Fe and As and a fractional Ba is not as 
straightforward as in the case of the 1111 material. Es- 
pecially the fact that two consecutive FeAs layers cannot 
be mapped by a mere translation in z direction leads to a 
dependence of part of the interlayer hopping parameters 
on the respective sublattice position, and eventually to 
the necessity of introducing additional kz dispersions. In 
the following 1 denotes the dxz^ 2 the dyz^ 3 the dx2_y2^ 
4 the dxy and 5 the d^z'^-r'^ orbital. The hopping param- 
eters are tabulated in Table |T] and Table [Til In addition 
we have the 4 onsite energies measured from the Fermi 
energy as ei/2 = 0.0987, 63 = -0.3595, 64 = 0.2078, and 
65 = -0.7516. 



^11/22 = 2t^/^ COS kx + 2t^/^ cos ky + 4^^ cos k^ cos ky ± 2t]^^{cos 2kx - cos 2ky) + ^t^y^^yy cos 2kx cos ky 

-^^tlly/xxy 2^2/ + "^tlxyy 2kx COS 2ky + 4:tl\{cOS kx + COS ky) COS kz 

±4t^^^(cos 2kx — COS 2ky) cos kz 
^33 = 2t^^(cos kx + COS ky) + 4t^^ cos kx cos ky + 2t^^(cos 2kx + cos 2ky) 

^44 = 2t^^(cos kx + cos ky) + At^y cos kx cos ky + 2t^^(cos 2kx + cos 2ky) + 4t^^^(cos 2kx cos ky + cos 2ky cos kx) 

'^^^xxyy '^^x COS 2ky + 2t^^ COS kz + 4t^^(cOS kx + COS ky) COS kz + ^t^yz COS kx COS ky COS kz 

^55 = 2t^^(cos kx + COS ky) + 2t^^(cos 2kx + cos 2ky) + 4t^^^(cos 2kx cos ky + cos 2ky cos kx) 
-^^t^xxyy COS 2kx COS 2ky + 2t^^ COS kz + 4t^^(cos kx + cos ky) cos kz 

<fi2 = 4t^^ sin kx sin ky + 4t^^^(sin 2kx sin ky + sin 2ky sin kx) + "^tx^xyy ^^^c sin 2ky + 8t^^^ sin kx sin ky cos kz 
^"13/23 = 2^4^ sin ky/x + 4:itll ^v/^ ^^/v ~ ^^^llyi^^^ ^K/x cos kx/y - cos 2kx/y sin ky/x) 

^14/24 = sin kx/y ± 4it^^ cos ky/x sin kx/y ± ^'^t]^xy si^ '^^x/y cos ky/x ± 4zt^^ sin kx/y cos kz — ^t^^z si^ ^x/y sin kz 

^^it]^y^ COS ky/x sin kx/y cos kz ± ^^txtyz sin 2kx/y cos ky/x cos kz — ^t^xyz sin 2kx/y cos ky/x sin kz 
^15/25 = ±2^4^ sin ky/x T 4^4^ sin ky/x cos kx/y T 8^4^^ sin ky/x cos kx/y cos kz 
^34 = 44^2/ (^^^ '^^y ~ 

^35 = 24^ (cos kx — cos ky) + 4i^^y{co^ 2kx cos ky — cos 2ky cos kx) 

^45 = 44^ sin kx sin ky + ^t'^xyy sin 2kx sin 2ky + 2^4^ sin kz + 4^4^ (cos kx + cos ky) sin kz 
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TABLE I. The intraorbital hopping parameters used for the DFT fit of the 5 orbital model. 



imm 




i = X 


i = y 


i = XX 


i — xy 


i — xxy 


i — xyy i — xxyy 


i — z 


i — xz i — xxz 


i — xyz 


m — 


1 


-0.0604 


-0.3005 


0.0253 


0.2388 


-0.0414 


-0.0237 0.0158 




-0.0101 0.0126 




m = 


3 


0.3378 




0.0011 


-0.0947 












m = 


4 


0.1965 




-0.0528 


0.1259 


-0.032 


0.0045 


0.1001 


0.0662 


0.0421 


m = 


5 


-0.0656 




0.0001 




0.01 


0.0047 


0.0563 


-0.0036 





TABLE IL The interorbital hopping parameters used for the DFT fit of the 5 orbital model. 



imn 




i = X 


i — xy 


i — xxy 


i — xxyy i — z i — xz 


i — xyz 


i — xxyz 


mn = 


12 




0.1934 


-0.0325 


0.0158 


-0.0168 




mn = 


13 


-0.4224 


0.0589 


0.0005 








mn = 


14 


0.1549 


-0.007 


-0.0055 


0.0524 


0.0349 


0.0018 


mn = 


15 


-0.0526 


-0.0862 






-0.0203 




mn = 


24 








0.0566 




0.0283 


mn = 


34 






-0.0108 








mn = 


35 


-0.2845 




0.0046 








mn = 


45 




-0.0475 




0.0004 -0.019 -0.0023 
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